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! Gemini North adaptive optics imaging spectroscopy is presented for the 

! Galactic ultra-compact H II (UGH II) region K3-50A. Data were obtained in 
the i^— band using the Near-infrared Integral Field Spectrograph (NIFS) behind 

^ ! the facility adaptive optics module ALTAIR in natural guide star mode. The 

r-| I NIFS data cube reveals a complex spatial morphology across the 0.1 pc scale of 

^1 the 3" UGH II region. Gomparison of the nebular emission to Gloudy ionization 

2 I models shows that the central source must have an effective temperature be- 

^ [ tween about 37000 K and 45000 K with preferred values near 40000 K. Evidence 

li^l is presented for sharp density variations in the nebula which are interpreted as 

I I a clearing of material nearest the central source. High excitation lines of Fe III 

' and Se IV show that the ionization of the nebula clearly changes with distance 

CO 

\ from the central source. A double lobed kinematic signature (± 25 km s~ ) is 

^ I evident in the Br7 line map which may be related to the larger scale ionized 

^ I flow detailed in earlier investigations. This signature and the large scale flow are 

O I not co-alighned, but they may still be related. Though clearly resolved from the 

^ I nebula, the central source itself remains buried, and the NIFS spectrum shows 

\ no evidence of photospheric lines. 
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1. Introduction 

K3-50A has been the subject of many studies in the optical/near-infrared and radio 
regimes and most wavelengths in between. It is highly obscured at short wavelengths and 
bright at long wavelengths. The object was originally part of the planetary nebulae survey 
of Kohoutek (1965), but was later reclassified as an H II region. Persson & Frogel (1974) 
observed K3-50A at mid-infrared wavelengths and showed the optical nebula was offset 
from the radio source due to a strong gradient in the foreground extinction. Wynn- Williams 
et al. (1977) determined precise offsets between the radio, infrared and optical sources which 
validate the Persson & Frogel (1974) model. 

K3-50A, also known as the Galactic radio source G70.3-1.6, has been extensively ob- 
served at long wavelengths. Turner &: Matthews (1984) provided a high angular resolution 
map at 2 cm which resolved the ionized gas into a clumpy arc-like or shell-like structure on 
scales smaller than 1". Kurtz et al. (1994) included this object in their catalog of UCH II 
regions, describing it as a core-halo type while De Free et al. (1994) observed the larger neb- 
ula in the continuum and radio recombination line H76q;. The radio morphology is extended 
in two lobes to the NW and SE and De Pree et al. (1994) found a velocity gradient along 
the lobes which led them to suggest the structure was an ionized outflow extending lO's of 
arcseconds from the center of the UCH II region. 

Howard et al. (1996) presented Br7 and Bra emission-line maps at greater than arcsec- 
ond resolution and followed this work (Howard et al. 1997) with molecular line maps which 
they modeled as arising from a large scale molecular torus (roughly perpendicular to the 
ionized outflow of De Pree et al. 1994). 

More recent work has centered on near-infrared and mid-infrared wavelengths and 
higher angular resolution data. Hofmann et al. (2004) observed the nebula in the if —band 
using speckle techniques to provide an order of magnitude higher angular resolution than 
previous near-infrared work, and Okamoto et al. (2003) used the Subaru 8m telescpe and 
mid-infrared imager COMICS to obtain ^0.5" resolution images and spectra near 10 //m. 

In this paper, we present if —band integral-fleld spectroscopy of K3-50A obtained on 
the Gemini North 8m telescope with the Near-infrared Integral Field Spectrograph. These 
data were taken in conjunction with the facility adaptive optics system and result in the 
most detailed data cube yet for the core of K3-50A. In an earlier paper (Blum & McGregor 
2008), we presented similar data for the UCH II region G45.45+0.6. In that case, several 
massive stars were identified by photospheric lines and other objects showed signatures of 
buried massive young stellar objects. G45.45+0.6 clearly harbors a rich cluster. Messineo et 
al. (2007) make similar conclusions based on their ESO VLT SINFONI integral-field spectro- 
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graph observations for the UCH II region in the cluster [DBS2003]8. Previous investigators 
have suggested the same may be true of K3-50A: Colgan et al. (1991), Okamoto et al. (2003), 
and Hofmann et al. (2004). Our data identify multiple continuum sources in K3-50A and 
support the multiple star model, but we also show that the excitation in the central few 
arcseconds of K3-50A is likely dominated by a single object (the unresolved continuum 
peak). 

In the following, the distance to the nebula is based on the work of Harris (1975). The 
Harris (1975) distance is derived from a compilation of velocities and a rotation model; 
Harris' distance needs to be reduced due to a revised distance to the Galactic center (Harris 
1975, used a distance of 10 kpc). Adopting a value of 8 kpc gives a reduction by a factor of 
0.8 resulting in a distance to the H II region of about 7 kpc. 

2. Observations 

Data were obtained with the Near-infrared Integral Field Spectrograph (NIFS) at the 
Cassegrain focus of the Gemini North Fredrick C. Gillett 8-m telescope on Mauna Kea, 
Hawaii on the night of 20 July, 2006 (HST). NIFS was used with the facility adaptive optics 
(AO) module ALTAIR ^ in natural guide star (NGS) mode. 

NIFS is fully described by McGregor et al. (2003); see also Blum & McGregor (2008) 
for a more detailed description of similar observations with NIFS as described here. Briefly, 
NIFS shces an approximately three arc second by three arc second field into 29 spectral 
segments of 0.1" width (the "slit" width) and ~ 3" in length. The scale along the spatial 
dimension is 0.043" pix~^. The resulting "spaxels," or spatial pixels, are thus 0.043" x 0.1" 
in size, and each contains a full spectrum covering one of the near-infrared bands. 

In the present paper, /T— band spectra are presented for the Galactic UCH II region K3- 
50A. The AO guide star used for ALTAIR is located 10.2" to the south and west of K3-50A 
and has an R magnitude of 12.7 according to the USNO catalog (source id: U1200_14213688). 
There was thin cirrus at the time of the observations, and the seeing (at 5000 A, corrected 
to zenith) reported during the observations by ALTAIR was approximately 0.45". The 
observations were obtained at relatively high airmass (1.5), and ALTAIR was run at 500 Hz. 

The observations consisted of a single coadded frame taken on source and a second frame 
obtained on a nearby (~ 50" West) blank field. Each frame had a total exposure time of 
600 seconds (15 x 40s) on source and 600 seconds on sky. NIFS was oriented with a position 



^http: / /www.gemini.edu/sciops / instruments / altair /altairlndex.html 
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(PA) of 70 degrees east of north to correspond approximately to the orientation of the large 
scale ionized flow from K3-50A (De Pree et al. 1994). This PA is aligned with the outflow 
direction. 

The spectral resolving power of NIFS in the band is A/ A A = 5160 which results in 
a linear dispersion of 2.13 A/pixel. This dispersion, combined with the large format array 
gives a full wavelength coverage at K of about 4200 A accounting for some minor truncation 
in the final data cube due to the systematic shift of wavelength in each slitlet from the 
staircase design of the image sheer. 

3. Data Reduction 

The line maps and spectra presented here were obtained using the Gemini NIFS IRAF^ 
data reduction package (version 1.9). 

The NIFS IRAF package allows for full reduction to the "image cube" stage where a 
final cube has a roughly 60x62 pixel image plane and a 2040 pixel spectral depth. First, raw 
images are prepared for reduction by standard Gemini procedures that create the FITS image 
variance and data quality extensions. Next, the data are sky subtracted, fiatfielded, rectified 
spatially, and wavelength calibrated (see Blum & McGregor 2008, for a full description of 
these steps). The last two steps result in a uniform spatial-spectral trace for each image 
slice row and a linear wavelength scale. It is important to remember that the final IRAF 
data cube resamples each shtlet to two 0.05" pixels for convenience; the angular resolution 
in this dimension is still 0.1". 

A lamp image of Ar and Xc lines was used to determine the dispersion along each 
slice and as a function of the spatial dimension of each slice. A quadratic polynomial was 
used in this case which produced typical uncertainties in the position of a spectral line of 
approximately ± 0.2 A (~ 1/10 pixel). Final spectra are interpolated to a linear wavelength 
solution. 

The spectra were next corrected for telluric absorption by division by the spectrum of 
an AO V star (HIP 102568). The spectrum of HIP 102568 was corrected for intrinsic Br7 
absorption by fitting a Voigt profile to the telluric spectrum. Particular care was taken in 
the fit in order to remove this feature from the telluric standard as accurately as possible. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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The final wavelength calibrated, spatially rectified, and telluric corrected images were 
transformed into data cubes. The spatial pixel scale was resampled to 0.05" in the fine 
dimension and block replicated to 0.05" in the course dimension providing for a uniform 
scale as described in the previous section. Line maps presented below were extracted from 
these cubes using NFMAP. 

The zero point of the wavelength calibration is confirmed by the positions of OH airglow 
lines in the sky image. Analysis of these lines shows the zero point is accurate to ± 2.5 km 
s~^. The calibration results in an observed blueshift of the nebular lines of approximately 
3.8 A (1.8 pixels or —53 km s~^). Lines arc identified on a spectrum extracted from a 2.5" 
diameter aperture centered on the NIFS FOV; see Figure L The spectrum is flux calibrated 
using the results of Howard et al. (1996, see below). The fine positions and relative strengths 
compared to Br7 for this integrated spectrum are given in Table 1. 



4. Results 

An image of the continuum near 21700 A is shown in Figure 2 along with the radio 
continuum image from De Free et al. (1994). The continuum passband extracted was 9 
pixels (19.2 A) wide. The appearance of sources in our continuum image is similar to the 
speckle band image presented by Hofmann et al. (2004) who report an angular resolution 
of 0.11" in their image derived from observations obtaddiained under 1.1" seeing. Hofmann 
et al. (2004) identify 10 point sources in their image, though some of these appear more 
extended than others, and this is consistent with the continuum image shown in Figure 2. 
The brightest source in Figure 2 is source ^^1 of Hofmann et al. (2004). The size of the most 
compact source in our image, source ^8 of Hofmann et al. (2004), is about 3.7 pixels in the 
fine sampling direction suggesting a FWHM of ^0.19". Source #1 is more extended than 
this in the present image (~ 4.9 pixels), but this is due to crowding with nearby blended 
sources and the intense nebular emission surrounding it. Hofmann et al. (2004) report Source 
#1 to be "point hke." 

The 10 sources identified by Hofmann et al. (2004) are clearly seen in Figure 2 with the 
possible exception of their source #2 which is close to source #1 and not obvious in our 
image. The Hofmann et al. (2004) sources are indicated by filled circles in Figure 2. An 
additional point source (named "8N") appears near —0.6, —1.2 which is indicated by a single 
contour in Figure 2 of Hofmann et al. (2004). In addition to these objects, there are other 
knots or clumps which might contain embedded sources, or may be density enhancements in 
the H II region. These clumps appear to match well the morphology in Figure 1 of Hofmann 
et al. (2004). An exception is the "X-like" structure centered on source #1 described by 
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Hofmann et al. (2004). We see no evidence of this feature in the lower contours of the nebular 
emission even though similarly bright clumps from the Hofmann et al. (2004) image do show 
up in our image. The mid-infrared images of Okamoto et al. (2003) also do not show this 
feature, though diffraction spikes from the telescope secondary are present in their images. 

Okamoto et al. (2003) identified four bright sources in the core region of K3-50A covered 
by Figure 2 in 9-13 //m images with somewhat lower angular resolution. The three brightest, 
OKYMl-3, were unresolved as individual point sources. 0KYM3 was associated with source 
#1 by Hofmann et al. (2004). The source labeled STHOl in Figure 2 is the peak of [Ne II] 
emission in the maps of Okamoto et al. (2003). These authors suggest STHOl is the location 
of an embedded massive star. There is no obvious point source in our continuum image or 
the band image of Hofmann et al. (2004) at this position. 

The diamonds in Figure 2 indicate the mid-infrared sources identified by Okamoto et 
al. (2003). The offsets between these sources (OKYM 1-4 and STHOl) are taken from 
Okamoto et al. (2003) and the postion of OKYM 3 is assumed to coincide with Source #1 
(also assumed to be the case by Hofmann et al. 2004). Hofmann et al. (2004) identified 
a band source with each of these Okamoto et al. (2003) sources within the positional 
uncertainties of the mid-infrared image. These sources are the obvious ones in Figure 2 near 
each of the OKYM sources. 



4.1. Line Emission Maps 

The continuum subtracted Br7 emission is shown in Figure 3. Howard et al. (1996) 
present narrow-band Brackett imaging. Their data are seeing limited and do not resolve 
any of the structure shown in Figure 3. Howard et al. (1996) report a Br7 fiux of 3.7x10"^^ 

erg cm~^ s^^ in a 3.4" aperture centered on the Br7 peak in their images. This covers 
approximately the NIFS field of view. Equating this fiux to the total counts in an approxi- 
mately 2.5" aperture results in a peak intensity (one pixel) of 5.5x10^^^ erg cm^^ s^^ pix^^ 
or 2.2x10"^^ erg cm~^ s~^ arcsec"^. The peak Br7 emission is to the SW of the continuum 
peak, although there is a second peak nearly as strong and very compact (point like, super- 
posed on extended emission) centered on the continuum peak (source #1). The Hofmann et 
al. (2004) sources #4-#7 are also centered on local Br7 peaks. The rest of the Hofmann et 
al. (2004) continuum sources are offset from nearby local Br7 emission peaks. 

He I emission follows the Br7 emission closely. The left panel of Figure 4 shows the 
ratio of He I to Br7 for the He I 3po_3g i^^^ 21127 A. The ratio is remarkably uniform 
over most of the FOV and has a value of approximately 0.04 (uncorrected for extinction). 
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This value, which would be larger in the absence of extinction, already suggests a hot star 
(s) with effective temperature ;^40000 K (Benjamin et al. 1999; Hanson et al. 2002). The 
stronger He I line, ^S-^P° at 20587 A, is not uniform, however, in comparison to Br7 (see 
the right panel of Figure 4 and discussion in § 5.2.3). 

Figure 1 shows that strong emission from multiply ionized species exist in the UCH II 
region. These include [Fe III], [Kr III], and [Se IV]. The first and last are the strongest 
lines apart from He I and Br7. Figure 5 shows the line maps for two of these lines. The 
[Fe III] emission is clumpy and peaks on the source 0KYM4 while the [Se IV] emission 
is more intense and somewhat more concentrated to the central "rectangle" or "diamond". 
The [Se IV] emission has a maximum at the continuum peak (—0.2, 0.6) which is unresolved 
(i.e. compact). The [Fe III] emission also has a local maximum at the continuum peak, but 
it appears to be affected by one or more bad pixels at the very center. These lines (including 
[Kr III]) appear to be common in UCH II; see Hanson et al. (2002) and Blum & McGregor 
(2008). The ionization potentials of Fe III, Fe IV, Kr III, Kr IV, Se III, and Se IV are 30.7, 
54.8, 37.0, 52.5, 30.8, and 42.9 eV, respectively. This suggests that the region with strong 
[Se IV] emission corresponds with Fe IV (for which there are no identified lines). The range 
of species seen is consistent with the presence of He I recombination lines (He I ionization 
potential = 24.4 eV) and the lack of He II recombination lines (He II ionization potential = 
54.4 eV). 

Molecular hydrogen emission is also detected (see Figure 1 and Table 1). The strongest 
line is the v=l-0 S(l) line at 21218 A. A map of the line is shown in Figure 6 with the 
contours of the Br7 emission overlaid. The projection of the H2 emission seems to follow 
the ionized emission with the peak emission coincident with the Br7 peak. There is a local 
maximum at the position of the continuum peak as well. The location and excitation of the 
H2 emission are discussed further in § 5.4. 

4.2. Gas Kinematics 

Figure 7 shows the measured Br7 line center velocity (in km s^^) for the entire FOV. The 
mean observed Br7 velocity of —50 km s~^ has been subtracted to show relative velocity 
differences in this map. This mean corresponds to a velocity with respect to the Local 
Standard of Rest (LSR) Vlsr = —62.5 ± 5 km s~^. The velocities are overplotted with the 
Br7 emission contours. Each spatial pixel in the Br7 map was fit with a Gaussian profile 
over a restricted velocity range of —250 to 50 km s~^. This choice was made so that the 
He I line at 21642 A (see Table 2) did not interfere with the continuum fit. The velocity 
map shown in Figure 7 was made with a modified version of the Euro3D program (Sanchez 
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2004) produced by one of us (PJM) and called NF3D. 

The Br7 velocity map shows a ~ ± 25 km s~^ range in velocity centroids. This should 
be compared with our FWHM velocity resolution of ~ 60 km s~^. Two Br7 emission clumps 
are apparent in Figure 7. These are separated by 1" — 2" and oriented approximately at a 
position angle of 45 degrees. A single Gaussian provides an adequate fit to the Br7 emission 
over the entire FOV. However, as we show in the following section, there is a broad blue 
wing to the emission near the continuum peak. In what follows, we call the single Gaussian 
fit to the bright component of the line "the principal component." When two Gaussians are 
fit to spaxcls near the continuum peak, the effect is to redshift the principal component by 
a further '-^ 5 km s^^ in this vicinity. 

The continuum peak is not located near any point of symmetry with respect to the 
velocity distribution. The nearest possible point source to the apparent center of the velocity 
distribution (Figure 2) is STHOl which was not detected in our continuum map. This source 
is only detected in the mid-infrared [Ne II] map of Okamoto et al. (2003) and not in any 
continuum image presented here or by Okamoto et al. (2003) or Hofmann et al. (2004). 

4.3. The He I Spectrum and a Broad Br7 Component 

The UGH II region K3-50A exhibits a rich He I spectrum in the A'— band. Table 2 
lists the He I lines detected and their properties. Using the on-line data base compiled by 
Peter van Hoof^; 13 lines are identified between 20000 and 22300 A. Other lines redward 
of 22300 A may be blended with the H I Pfund series (Figure 1). Multiple He I lines are 
present in the region around B7. Figure 8 shows spectra extracted from individual spaxels 
at a typical location ~ 1.3" SE of the continuum peak (Figure 1, left panel) and a location ~ 
0.2" SE of the continuum peak (Figure 1, right panel). Five Gaussian functions have been fit 
simultaneously to these features using the spectral analysis program Liner^. Good fits with 
RMS residuals of < 1% have typically been achieved by assuming that all five lines have the 
same FWHM as Br7 (Figure 1, left panel). 

Poorer fits with RMS residuals of ~ 3% are obtained in the region within ~ 0.5" of the 
continuum peak where the Br7 line has a strong blue wing (e.g., Figure 1, right panel). The 
Gaussian component fit to this wing has a central wavelength (which is a fit parameter) that 



^http://www.pa. uky.edu/~peter/atomic/; the on-line database uses data from Ralchenko et al. (2007). 

^Liner is provided by R. Pogge at the Ohio State University. We appreciate Dr. Pogge's support of the 
Liner package. 
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is close to that of the He I 4/ - 7d ^F°-^D, ^F°-^D 21655 A hne (see Table 2). However, 
this line is expected to be too weak to account for the blue wing near the continuum peak. 
In fact, based on atomic data used in the Cloudy photoionization code (Ferland et al. 1998), 
the 21655 A He I hne is expected to be weaker than the He I 4/ - 7g ^F°-^G, '^F°-^G 21647 
A line that is apparent in Figure 1. Better fits are obtained for spaxels near the continuum 
peak if the FWHM of the blue wing component is allowed to vary. RMS residuals of < 1% 
are then obtained with a Gaussian FWHM that is about twice as broad as the principal 
Br7 hne. There is no apparent reason why the 21655 A He I hne should be anomalously 
strong and have a width twice that of other He I hnes. We therefore conclude that the blue 
wing on Br7 is due to a mix of a second kinematic component of Br7 that dominates in a 
small region near the continuum peak and weaker He I emission that occurs over most of the 
FOV. It is interesting that neither the peak emission nor the width of this second kinematic 
component is centered on the continuum peak. Instead, this emission forms an arc that is 
offset SE of the continuum peak by < 0.5". 

Whatever the source of the broad component, the contribution of the lines near 21600 
A make a significant contribution to the integrated "Br7" flux. For the spaxel shown in 
Figure 8, the broad line produces 50% of the corresponding Br7 flux. 

5. Discussion 

5.1. Source Positions &; Morphology 

Wynn- Williams et al. (1977) measured the optical, near-infrared, and mid-infrared 
positions of K3-50A with respect to the radio continuum (CoUey & Scott 1977) using 
background stars as a reference. The CoUey & Scott (1977) position is consistent with later 
radio positions including the high angular resolution study of Turner & Matthews (1984), and 
thus, to our knowledge, the Wynn-WiUiams et al. (1977) positions are still the best available. 
More recent studies have assumed the near and mid-infrared peaks are coincident with the 
radio continuum peak (Okamoto et al. 2003; Hofmann et al. 2004). Wynn- Williams et al. 
(1977) flnd that the 2.2 /im and 10 peaks are coincident, within the ~ 1" uncertainties, 
with the radio continuum peak. The optical and 1.6 /xm peaks are located to the south 
west of the radio peak, and this is consistent with an extinction gradient across K3-50A as 
suggested by Frogel & Persson (1974) and Wynn-WiUiams et al. (1977) and found also by 
Howard et al. (1996). 

Dc Free et al. (1994) observed K3-50A in the radio continuum (15 GHz) and in the 
H76a recombination line. The angular resolution of these data is about 1.5". The continuum 
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appears to be extended along an axis with PA (E from N) of 160 degrees and there is a strong 
velocity gradient from NW to SE along the extended continuum source (see also Roelfsema et 
al. 1988). This has led to the accepted model of an ionized outflow along the 160 degree axis. 
The NIFS continuum is shown superposed with the De Pree et al. (1994) radio continuum 
image in Figure 2. The two images have been aligned assuming the radio and iC— band 
continuum peaks are coincident. The large difference in scale and angular resolution are 
apparent. The overlay is instructive for the following discussion, but there are no obvious 
connections which can be drawn between features in the two maps. 

Molecular line studies have been used to argue that there is a rotating torus of material 
with major axis roughly coincident with and perpendicular to the radio source. Low angular 
resolution CO data and modeling of HCO+ line emission suggests this material may be 
rotating about an axis parallel to the outflow axis. De Pree et al. (1994) first drew attention 
to a possible link between their ionized outfiow and the molecular gas traced by CO line 
emission (Phillips & Mampaso 1991), though they cautioned that the low angular resolution 
of the CO data (30") and the marginal evidence for a velocity gradient along its major axis 
made any link between the two structures "premature." Later, Howard et al. (1997) used 
observations of HCO^ to argue for a rotating torus roughly perpendicular to the ionized 
outflow. The peak lobes of HCO"*" emission are separated by 8", so the scale of the rotating 
molecular torus is larger than the FOV covered by NIFS. The velocity difference measured 
between the peaks of the HCO+ is ~ ± 2 km s~^, and the best fit to the velocity profile was 
found to be proportional to the square root of the radius. 

Higher angular resolution observations were made by Okamoto et al. (2003) in the 
mid-infrared (near 10 /xm) and by Hofmann et al. (2004) in the band (2.2 /xm). The 
former observations were diffraction hmited based on the 0.4" FWHM obtained on standard 
stars, but the images of K3-50A do not clearly separate individual objects owing to intense 
resolved structure/emission and blended point sources. Okamoto et al. (2003) fit various 
point sources (Figure 2) to their data. The Hofmann et al. (2004) images, on the other hand, 
were obtained with a speckle imaging technique and are of similar angular resolution as the 
data presented here. The higher angular resolution and lower background emission result in 
more identifiable point sources within the central few arcseconds of K3-50A. Okamoto et al. 
(2003) and Hofmann et al. (2004) conclude that K3-50A contains multiple massive stars; the 
former through excitation measurements of mid-infrared ionized lines and the latter based 
on the detection of multiple point sources in the center of the UCH II region. 

The NIFS data confirm various aspects of the earlier work, principally the /T— band 
morphology of the point sources and brighter clumps of emission which are dominated by 
Br7. The band image of Hofmann et al. (2004), the Okamoto et al. (2003) mid-infrared 
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images, and our NIFS images all indicate an abrupt cutoff in emission to the north of the 
bright continuum peak. The gradient is extremely sharp and the emission drops precipitously 
in ^0.5" towards the north of the UCH II region. The gradient in extinction was noted in 
earlier work (e.g., Wynn- Williams et al. 1977) but is seen very clearly in the J, if, and K 
images of Okamoto et al. (2003). The center of the emission shifts systematically to the 
south. This geometry has been described as an opening cone with the continuum peak at 
the cone's apex (Hofmann et al. 2004) and the opening toward the south-south east. 

The dense molecular torus (Howard et al. 1997) and other overlying material form the 
cone structure by obscuring ionized material to the north; sec Figure 7 of Howard ct al. 
(1997). The constraining effect of the torus may be responsible for helping to shape the 
ionized outflow (De Free et al. 1994). The very strong gradient in emission evident in 
Figure 2 and the band image of Hofmann et al. (2004) suggest a possible obscuring disk 
on small scales, but no kinematic or morphological evidence exists to support this picture, 
and the radio and near-infrared peaks are offset from the center of the torus (see Figure 14 
of Okamoto et al. 2003). If the embedded massive star associated with the continuum peak 
formed from the molecular torus, then it may have separated itself from this "parental cloud." 



5.2. Ionization of the Nebula 

Okamoto et al. (2003) presented ionic line measurements of coUisionally excited Ar 111, 
S IV, and Ne II at approximately 9, 11, and 13 //m, respectively. Ratios of these lines suggest 
a relatively low excitation, typical of late 0-type or early B-type stars. The radio continuum 
(Turner & Matthews 1984; De Free et al. 1994), however, requires more ionizing photons 
than a single late 0-type star produces. Okamoto et al. (2003) thus concluded that multiple 
stars ionize K3-50A (as discussed above for the NIFS images, we detect continuum sources, 
#1, #4-#7 at least, which are centered on local peaks of the Br7 emission). Okamoto 
et al. (2001) made a similar analysis and conclusion for W51 IRS2, also a massive star 
forming region. But recently, Barbosa et al. (2008) identified the exciting star of W51 IRS2 
directly through its near-infrared spectrum, and it is much hotter (approximately 03-04). 
For G29.26— 0.02, Morisset et al. (2002) find a best fit nebular model which also has a 
cooler central source than the observed spectral type (Watson & Hanson 1997), though their 
acceptable models do cover a range of Tcfj which includes the observed spectral type. The 
nebular emission in G45.45+0.06 was consistent with the observed spectral types found by 
Blum & McGregor (2008), but detailed modeling was not done. 

A grid of ionization models has been run to analyze the near-infrared lines presented 
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here in order to further explore the excitation of the UCH II region. Cloudy^ was run on a 
grid of input parameters which sample a range of ionizing flux and T^s of the central source 
as described by Blum & McGregor (2008). A base line model grid was run with density 
10^ cm~^ and the Castelli & Kurucz (2004) atmospheres (log (g) = 5.0 and Z = 0). The 
simulations were run until the temperature into the cloud reached 30 K. A model He atom 
was used with 10 energy levels to better resolve the level populations for the near-infrared 
hues. The " Jsm" abundance set available in Cloudy was used which has a He abundance by 
number relative to H of 0.10. In the following sections, we discuss how the NIFS observations 
constrain inputs to the models and also what can be inferred from a comparison of the models 
and observations. 



5.2.1. [Felll] Lines 

The density used in the Cloudy runs was constrained by the observed ratio of [Fe III] 
hues. The ratio of [Fe III] 22184 A to 21457 A is sensitive to density (Lutz et al. 1993; 
Bautista & Pradhan 1998). The observed ratio (see Table 1) is 3.12 and corrected for A^ 
of 2.17 (Howard et al. 1996) leads to a ratio of 2.8 and a density of approximately 10^ cm~^ 
according to the 7000 K model of Lutz ct al. (1993) or the 10000 K model of Bautista & 
Pradhan (1998). Bautista & Pradhan (1998) give densities for the ratio of [Fc III] 22427 A to 
22184 A as well, but the observed ratio of 0.5 is on the fiat part of their curve and gives a low 
density of ^10^ cm~^. The small reddening correction would push toward even lower density. 
All three lines are stronger toward the south in the FOV and weak near the continuum peak, 
thus the ratios do not constrain the density near the continuum peak. Turner & Matthews 
(1984) derived a density of 10^ cm"^ in the central few arcseconds from their 2 cm data; 
similar radio data were presented by Kurtz et al. (1994) who find a density of 4x10^ cm~^. 
We discuss models for a range of densities below (see §5.2.2 and §5.2.3 and Figure 11). 

The [Fe III] line ratios to Br7 reported in Table 1 are similar to those reported by DePoy 
& Pogge (1994) in Orion. Though the emission is relatively strong in the K3-50A spectrum, 
it may be consistent with most of the Fe being condensed onto dust grains as is the case 
for Orion; similar ratios and physical conditions in Orion are due to a gas phase abundance 
which is only a percent of the solar iron abundance (DePoy & Pogge 1994). The distribution 
of Se IV and Fe III discussed above suggests some of the gas phase iron in K3-50A must be 
in Fe IV, but there are no Fe IV lines to estimate the abundance in that species (the upper 



^Calculations were performed with version 07.02 of Cloudy, which was last described in detail by Ferland 
et al. (1998). 



-13- 



levels of observable lines are generally too high to be excited in H II regions (Bautista & 
Pradhan 1998)). For the models described below which fit the data best for the excitation of 
the nebula, most of the Fe is in Fc III and Fc IV. In some cases, the Fe IV fraction is ~ 90 %, 
but this still requires the bulk of the iron to be condensed on grains if the total abundance 
is near or above the solar value. 



5.2.2. Recombination Lines and Excitation 

The map of the observed ratio of He I 21127 A to Br7 is shown in the left panel of 
Figure 4. The ratio is approximately uniform across the main emission in the center of 
the image and has a value of between 0.04 and 0.05. The highest value is 0.055 near the 
continuum peak, and regions around the periphery have values in the range 0.025-0.03. The 
value in the average spectrum of Figure 1 is 0.046 (Table 2) which includes the effect of the 
blend of lines near Br7. These other lines are not included in the model prediction shown 
in Figure 9, although their effect is small. The average ratio of 0.046 is consistent with an 
0-star which has Tefr;^40000 K (Hanson et al. 2002). A small correction for extinction gives 
an intrinsic ratio of 0.051 if we adopt the observed ratio in Table 2. 

Cloudy models were run for a grid of ionizing flux (phi, s~^ cm^) versus exciting source 
Teff . The former parameter was varied from 12.0 to 16.0 dex and the latter varied from 25000 
K to 45000 K. The results for our basehne case (density = 10^ cm~^ and bhster geometry) are 
given in Figure 9. The ratio along a given contour becomes approximately constant above 
Teff = 40000 K, but the detailed ratio depends on the ionizing flux as well. In Figure 9, the 
overall emission due to dust compared to the Br7 flux as predicted by Cloudy is also plotted 
(as dashed contours). The intersection of this contour with the He I to Br7 contour can give 
an idea of the true excitation of the nebula under the assumption of a single dominant source 
(which might not be the case); see Blum & Damineli (1999) and Blum & McGregor (2008). 

The dust emission was calculated approximately from the spectral energy distribution 
(SED) shown in Figure 2 of Howard et al. (1997). No attempt was made to correct for 
the different beam sizes associated with the many distinct data sets which make up the 
Howard et al. (1997) flgure. The most important is the ~ 100 /xm point originally reported 
by Thronson & Harper (1979) since this is the peak in the SED. The Thronson & Harper 
(1979) beam was 50" in diameter and included flux from the nearby K3-50B H II region. On 
the other hand, higher angular resolution data at 450 /im (Thompson et al. 2006) show a very 
compact source (unresolved at the 8" beam size) with a similar flux to the SED presented 
by Thronson & Harper (1979). Howard et al. (1996) report a Br7 flux of 3.2x10^^^ ergs 
cm~^ s~^ corrected for Ak = 2.17 mag. We flnd a integrated dust emission of approximately 
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1.8x10 ^ ergs cm ^ s ^ from the Howard et al. (1997) SED. Thus the lower hmit to the 
ratio is 1.78x10-^ 

The intersection of the contours for He I and dust compared to Br7 suggest a T^s of 
about 38000 K and LoglO(phi) = 13.8. Since the Br7 to dust ratio is a lower limit according 
the discussion above, the true Tes would be hotter. The ionization model grid was repeated 
for a density of 10^ cm~^. The results are similar to those in Figure 9 with the intersection of 
the contours for the observed ratios moving to slightly lower T^s and higher ionizing flux, i.e 
up and to the left in the figure. The formal intersection is at 37000 K and Logl0(phi) = 14.5. 
Changing the geometry from plane parallel (open) to spherical (closed) in Cloudy requires 
higher incident ionizing flux due to the radiation field dilution by In the model grid, 

the inner radius of the cloud was set to 0.5" or 5x10^^ cm. The observed ratios for the 
spherical case intersect at somewhat higher Teflf, 40000 K, and the corresponding ionizing 
flux is LoglO(phi) = 14.7 to 14.4. If the dust emission is too high by a factor of two due to 
beam size mismatches, then the Br7 to dust ratio is too low by a factor of two. In this case, 
the model Teflf could be as high as 45000 K for a spherical geometry. 

We conclude that the excitation of K3-50A is consistent with a central star whose Teflf is 
greater than 37000 K and less than 45000 K. This is higher than Okamoto et al. (2003) derive 
based on line ratios of collisionally excited ions. A similar result is obtained in W51 IRS2 
when comparing the results of Okamoto et al. (2001) for line ratios to the recent spectral 
type obtained for the star located at the position of 1RS2 (whose UCH II region is known as 
W51d); see Barbosa et al. (2008). Evidence is presented below that the geometry may be 
more open like, see also] []col91, depree94 which would lead to Tcs values near 40000 K. 

De Free et al. (1994) report a Lyman continuum luminosity, LyC = 2.1x10^^ s~^ for a 
distance of 8.7 kpc. Adopting 7 kpc (see §1), results in a reduced number of 1.3x10^^ s~^. 
Kurtz et al. (1994) find LyC = 1.95x10^^ s~^ for a distance of 8.3 kpc, which scaled to 7 kpc 
gives a similar number. This LyC luminosity is consistent with a single 0-type star with 
effective temperature of approximately 40000 K based on the properties for O stars given 
by Martins et al. (2005) or 41000 K if the data from Vacca et al. (1996) are used. Figure 2 
shows there are multiple continuum sources (especially in the vicinity of the peak emission) , 
but the excitation of the nebula based on the NIFS spectrum is consistent with a single, 
dominant star. 
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5. 2. 3. He I 20587 A Emission 

The strongest He I line in K3-50A is the 20587 A hne (Figure 1). The right panel 
of Figure 4 shows the line ratio relative to Br7. Comparison to the left panel of Figure 4 
indicates the 20587 A line emission is more variable across the nebula than is the case for 
the 21127 A emission. The strong degree of uniformity for the latter suggests the 20587 
A line variations are not due to differential extinction, but instead are due to structure in 
the nebula. If the variations seen in the right panel of Figure 4 were due only to variations 
in the line of sight extinction, then this would require a change in Ak of more than 10 
magnitudes. The continuum map morphology is not consistent with this picture, nor is 
the average extinction found by Howard et al. (1996) in the central few arcseconds (A^ ~ 
2). Finally, the Br7 and He I 21127 A intensities are near maxima at the position of the 
continuum peak while the 20587 A line itself is weak there. 

The formation of the He I 20587 A line is complicated and depends on the physical 
conditions in the ionized gas (Shields 1993; Ferland 1999)*^. A large ratio to Br7 is due 
primarily to high density (Shields 1993) which allows for enhanced collisional population 
of the 2 level. However, this level population can also depend strongly on He I La 
584 A resonance fluorescence which is affected by dust absorption of the He La hne and 
ionization of neutral H by the He La line since both processes destroy the photons which 
would otherwise pump the line (Shields 1993). The biggest changes in 20587 A line formation 
are expected for a range of central star Teg. The Shields (1993) models were for an average 
over the nebula; in the present case, the nebula is spatially resolved and clearly is dominated 
by a point source (or unresolved sources) at one end of the FOV. In Figure 10, the ratio 
of He I to Br7 emissivity with depth into the cloud (for a blister model) is plotted for the 
21127 and 20587 A lines and a central source with Teff = 38000 K. The 20587 A hne strength 
changes sharply from inner (close to central source) to outer (far from the central source) 
positions in the cloud. 

Thus, at any given point in the nebula, we would expect a relatively high ratio in the 
right panel of Figure 4 since the effects shown in Figure 10 would be integrated along the 
hne of sight. The range of observed values for the ratio in the right panel of Figure 4 then 
suggests that real density variations may exist across the field of view (lower density near 
the continuum peak) . A possible exception to this would be if we were looking along a line 
of sight such that a slab or bhster was illuminated with a geometry whose normal to the slab 
was in the plane of the sky. But this geometry is inconsistent with the molecular hydrogen 



^Thc numerical values of the 20587 A line have changed since the Shields (1993) paper; see Blum & 
Damineh (1999); Ferland (1999) 



-16- 



morphology (see below, §5.4) which suggests a face-on bhster geometry. 

In Figure 11, the cloudy predictions for the ratios of the 20587 A and 21127 A lines 
to Br7 versus density from Cloudy are plotted. The latter line has an observed ratio of 
approximately 0.051 (corrected for extinction) and is nearly uniform over the nebula. The 
former line has an observed ratio (corrected for extinction) that varies from about 0.33 to 
0.84. Figure 11 shows the Cloudy prediction for the case of T^g = 38000 K and LoglO(phi) 
= 14.0. These results indicate that the 21127 A line ratio is nearly constant over a range 
in density of 100 to 100000 cm^'^ while the 20587 A line changes from 0.38 to nearly 1.0 
over the same range. The ratios plotted in Figure 11 are integrated along the line of sight 
through the nebula in contrast to the ratios plotted in Figure 10. 

5.3. Kinematics 

The kinematics evident in the bright Bv-f line (Figure 7) do not fit neatly with any of 
the previous kinematic data that were obtained at lower angular resolution. Figure 7 shows 
a double lobed structure, possibly a fiow, with a major axis distinct from any of the larger 
scale components discussed above. Furthermore, the ordered red and blue shifted material 
is not well correlated with much of the strong Br7 emission and the point of symmetry of 
the indicated flow is not coincident with the continuum peak. A possible explanation for 
the small scale flow is that it is created by a lower mass protostar in the nascent cluster. 
This molecular outflow is then ionized by the massive star at the continuum peak. A similar 
scenario was recently described for the compact high mass star forming region IRS2 in W51 
(Lacy et al. 2007; Barbosa et al. 2008). 

Another possibihty is that the large scale ionized flow is related to the inner flow. This 
would require that the outflow is precessing or that details of the interaction of the out- 
flowing gas on larger scales than the NIFS FOV reorient the flow since the axis of the inner 
flow does not line up with the axis of the large scale ionized flow. The contours of the larger 
flow do, however, appear to twist (De Free et al. 1994) and Figure 2 clearly shows more 
detailed radio (kinematic) observations would be helpful in understanding the transition 
between the inner scale and larger scale gas kinematics. De Free et al. (1994) argued that 
the large scale ionized flow might result from the expansion of material coming from within 
a small (few arcsecond) shell. They took as evidence for a shell the high angular resolution 
2 cm map of Turner & Matthews (1984). The NIFS Br7 morphology appears to match well 
with the 2 cm emission, but does not appear shell like (Turner & Matthews 1984, noted 
the classiflcation of K3-50A as shell like was problematic). Still, the NIFS data suggest 
an outflow originating within this small region which is likely confined or directed by the 
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material around it. The large scale flow has a velocity gradient of about 6 km s"^ arcsec"^ 
while the small scale flow gradient evident in Figure 7 is ;^25 km s~^ arcsec"^. The velocity 
reported for the core of the radio emission is —33 km s~^ with respect to the LSR (De Free 
et al. 1994). The integrated Br7 emission has Vi^sR = —62.5 km s^^ (sec Table 1). The 
red-shifted emission in Figure 7 is at Vlsr ~ —30 km s"^. This suggests that we may be 
only seeing part of radio-emitting gas in the Br7 map and that the red lobe of the Br7 
flow is really tracing the bulk of the radio flow; that velocity is presumably the systemic 
velocity of the radio UCH 11 region. In this model, all the rest of the Br7 emission is tracing 
blue-shifted gas outflowing from the central source. The spatially unresolved Br7 emission 
at the continuum peak (see Figure 3 and §5.6) also has a Vlsr ~ —32 km that is close 
to the systemic velocity deflned by the H76q; emission. 



5.4. Nebular Structure 

In §4.3, a broad component of emission was identifled to the blue of the Br7 line (see 
Figure 8). The position of the line is consistent with one of the He 1 lines from the n= 7-4 
complex, but the emission is much stronger than expected based on transition probabilities 
and broader than any other He 1 line by about a factor of two. Thus, it is more likely 
a second Br7 component. The line velocity is about 40 km s~^blue shifted compared to 
the background or cloud velocity (the green regions in Figure 7). This is larger than the 
blueshifted velocity structure of the principal component of Br7 shown in that Figure. The 
width of the broad component is about twice the width of the principal Br7 component (120 
km s~^ compared to 70 km s~^). The broad component is distributed in an arc to the SSE 
of the continuum source (~ 0.5" in length and offset from the peak by ~ 0.2"-0.3") and is 
conflned to this vicinity of the continuum source. 

It is not clear how the high velocity, broad component is produced. It may be related 
to hot gas escaping a tightly constrained volume near the central source. The scale associ- 
ated with the broad component is small, only 0.007-0.010 pc (for a distance of 7 kpc). In 
champagne flows that are produced as the ionizing radiation from the hot star penetrates 
into a uniform medium (e.g., Yorke et al. 1984), the hot star quickly carves out a larger 
diameter cavity (;^0.1 pc). Perhaps shocked gas initially flowing parallel to the plane of the 
sky and then being redirected along our line of sight could explain the broad component. 
This putative flow might be related to the central source stellar wind. 

In any case, gas must escape to the north and eventually form the larger ionized flow 
seen in the De Free et al. (1994) radio images. The immediate vicinity of the continuum 
peak may be flUed with low density gas, and this gas may surround (or overlay) higher 
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density material within the projected "cone" which extends to the south of the continuum 
source. In order to explain the low ratio of He I 20587 A to Br7 near the continuum peak, 
the particular line of sight must contain very little dense ionized gas indicating this area has 
been cleared away as the hot star emerges from its parent cloud (which remains mostly to 
the north). The high angular resolution 2 cm image of Turner & Matthews (1984) and the 
molecular maps of Howard et al. (1997) are consistent with clearing of high density material 
in the vicinity of the continuum peak; both show voids or gaps in the emission. None of the 
current absolute positions is good enough to place the different wavelength data to better 
than about 0.5" which is the scale of the structure in our data and the 2 cm map. 

The nebula is seen in projection against H2 emission, primarily in the lower left quadrant 
where the Br7 emission is strongest (Figure 6) but also (unresolved) at the position of the 
continuum peak. Three prominent lines of molecular hydrogen are listed in Table 1. The 
average ratios of the two longer wavelength hues to the v=l-0 S(l) 21218 A hne (0.42 ± 
0.1, 0.23 ± 0.1 for H2 v= 1-0 S(0) and H2 v= 2-1 S(l), respectively) can be compared to 
the diagnostic diagram given by Hanson et al. (2002) which shows these same ratios for a 
large sample of UCH II regions along with several model predictions for shock excited and 
fluorescent emission. The bulk of the Hanson et al. (2002) sample is consistent with dense 
(10^ cm~^) photo-dissociation regions (PDR). The results for K3-50A place it in the same 
region of the diagram showing it too is consistent with a dense PDR on the boundary of its 
ionized nebula. Since we see both the ionized gas and molecular emission, it is likely the H2 
emission is behind the ionized gas (a large column of molecular gas in front would obscure 
our view). If this is the case, then the UCH II region may be best described as bhster like; 
the hot star has emerged from the parent cloud and is ionizing the face of the molecular 
material to the south. In this case, the "cone" is really a flattened structure seen face on 
and lying on the surface of a dense cloud of molecular gas. 



5.5. The Ne, Ar, and S Lines 

As mentioned above, Okamoto et al. (2003) analyzed [Ar IIl]/[Ne II] and [S IV] /[Ne II] 
ratios in K3-50A. The observed ratios are consistent with other UCH II regions, but in- 
consistent with simple ionization models. The models predict stronger [S IV] /[Ne II] than 
[Ar IIl]/[Ne II] for a given T^s, but the observations show the average ratios are of similar 
value and the peak ratios have substantially larger [Ar IIl]/[Nc II] than [S IV] /[Ne II]. Our 
grid of models agrees with Okamoto et al. (2003) in the sense that for a given T^s, the [S IV] 
line ratio is predicted to be stronger. 
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Martm-Hernandez et al. (2002) computed Ar, S, and Ne abundances for K3-50A from 
ISO mid-infrared spectra and find abundances which differ from the solar values (and our 
Cloudy models). In particular, the Ne abundance compared to H is approximately solar, 
while the Ar abundance is a factor of 2 lower than the solar value, and the S abundance is 
a factor of 10 lower. However, this large reduction in S abundance may be (in part) due 
to the high density in K3-50A depopulating the upper levels producing the S emission as 
Martin- Hernandez et al. (2002) point out. Detailed models used to derive the S abundance 
in K3-50A by Afflerbach et al. (1997) give a value only a factor of two lower than solar. The 
ISO spectra were obtained through large apertures (~ 20"), and the line ratios are smaller 
than the peak ratios in the higher angular resolution Okamoto et al. (2003) data. The ISO 
line ratios are nearly equal and have a value of approximately 0.25. This is comparable to the 
low end of the range of values obtained by Okamoto et al. (2003). The observed high peak 
hue ratios of Okamoto et al. (2003) ([S IV]/[Ne II] = 3, [Ar IIl]/[Ne II] = 4) would appear 
to require lower S abundance and higher Tcs to be consistent with the Cloudy ionization 
models, but the precise line ratios do not converge in our grid of models even when adopting 
the detailed abundances of Martm-Hernandez et al. (2002). In particular, reproducing the 
high hne ratio of [Ar III]/[Ne II] requires to be 45000 K, the hmit of the range based 
on the analysis of the He I 21127 A line and the dust emission, and such a model would 
produce a much too high [S IV] /[Ne II] ratio. 

Earher work based on Kuiper Airborne Observatory observations (Colgan et al. 1991) 
also did not succeed in fitting the ratios of these infrared lines. The Colgan et al. (1991) 
data were obtained through apertures approximately 5" in diameter and the line ratios are 
larger than the ISO values (about unity for both) and somewhat smaller than the Okamoto 
et al. (2003) peak values. Colgan et al. (1991) and Okamoto et al. (2003) have argued that 
multiple ionizing sources might explain the observations. 

The NIFS data generally (continuum. He I lines, resolved nebular structure) seem to 
indicate the central few arcseconds are dominated by the continuum peak. Figure 5 shows 
that the excitation in the nebula is higher toward the continuum peak. The [Fe III] emission 
is weak in the vicinity of the continuum peak compared to the [Se IV] line. The [Kr III] 
line is weak as well and this is consistent with its lower ionization potential than the other 
species present. In the case of 0KYM4 (Figure 2), there may be some evidence for lower 
excitation; Okamoto et al. (2003) preferred a late O or B star for this source. The left panel 
of Figure 4 shows a dip in the He I 21127 A to Br7 ratio at the position of 0KYM4 (the 
ratio is about 0.025). This source is about 2" away from the continuum peak and this is 
roughly the scale of an associated Stromgren sphere for the central star (Figure 10). Further, 
the [Fe III] emission peaks at the location of this source. Indeed, the emission morphology 
of the line maps suggest the boundary of the H II region which is principally due to the 
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continuum peak is on this same scale. If this is correct the "cone" morphology arises from 
dust obscuration to the north and the ionization boundary of the nebula to the south. The 
main nebula is largely excited by the continuum peak, and 0KYM4 may be locally internally 
excited. The compactness of 0KYM4 then suggests the source is very young and/or highly 
constrained. 

Lower mass OB stars should not greatly affect the structure within the "cone," however. 
Colgan et al. (1991) and Okamoto ct al. (2003) argued that the low average [Ar IIl]/[Ne II] 
and [S IV] /[Ne II] ratios might be due to multiple lower T^s sources (as may be the case for 
0KYM4), but a successful model must also explain the detailed ratios including their very 
high peak ratios near the continuum peak. While we have not been able to produce such a 
model, higher ratios are more consistent with hotter stars. 

5.6. The Continuum Pecik; 

In the preceding sections, the structure of the nebula was described, and it was concluded 
that the central source is a hot star emerging from its natal cocoon. The star is ionizing the 
cloud to the south in our FOV and produces an intense UGH II region and associated PDR. 
It was argued that the density in the nebula may be lower in the vicinity of the continuum 
peak. 

The continuum peak (by definition) dominates the i^— band emission. In addition, Br7, 
He I 21127 A, [Se IV] (the highest ionization hne observed), and H2 show unresolved emission 
peaks (or near peaks) at this location. There is clearly a very compact source of emission 
here which includes stellar continuum (from multiple sources) re-radiated by hot dust, and 
possibly line emission from the circumstellar environment. An ionized circumstellar disc 
is a possible source of the compact emission. Though no kinematic evidence is available 
to support this possibility. Bik & Thi (2004) and Blum et al. (2004) presented kinematic 
evidence of discs around (somewhat less) massive young stars based on the emission profiles 
of the CO v=2-0 rotational-vibrational bandhead at 22935 A. This line is not detected in 
K3-50A. Whatever its geometry, there must be compact circumstellar material (partially) 
surrounding the central source which is not directly linked (at least now) to the larger 
resolved UCH II region covered by the NIFS FOV. 
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6. Summciry 

High angular resolution (~ 0.2") i^'— band spectral imaging has been presented for 
the Galactic UGH II region K3-50A. The data were obtained with adaptive optics and 
the image slicing integral-field spectrograph, NIFS, at Gemini North. The combination 
of excellent wavelength coverage, moderate spectral resolution (5160), and high angular 
resolution provide the most detailed look to date at this massive young stellar object and 
its immediate environment. The emission morphology matches well with earlier speckle 
observations. 

Spatial and ionization structure are resolved in the 3" nebula. The He I emission 
morphology suggests changes in density across the field which are attributed to the break 
out and clearing phase of massive star birth. The ratio of He I 20587 A to Br7 varies from 
about 0.3 to 0.8 and requires that the density near the central source be more than a factor 
10 less than farther away. The ionization in the nebula can be traced across the nebula by 
the varying intensity of [Fe III] and [Se IV] lines. 

A grid of ionization models has been run with the code Cloudy (Ferland et al. 1998), 
and the results of these calculations were compared to the observed He I 21127 A line to Br7 
line ratio and dust emission properties. The combination of observed values constrain the 
nebular excitation to a source with T^s ^ 37000 K. This is somewhat hotter than expected 
from coUisionally excited emission-line ratios in the mid-infrared (Okamoto et al. 2003), but 
those line ratios, particularly their peak values on scales of ~ 1", are not explained by our 
models (or those of Okamoto et al. 2003). 

Despite the high angular resolution of these data which enhances the contrast between 
background and point source emission, no photospheric lines are detected in the spectrum 
of the bright continuum source; it remains buried within intense continuum emission. The 
central source is likely also blended with other point sources as pointed out by other in- 
vestigators. There is spatially unresolved, bright emission in the lines of Br7, He I, and 
[Se IV] (and possibly [Fe III]) at the continuum peak position and these are attributed to 
the circumstellar environment of the massive young stellar object. 

The NIFS data cube reveals a striking kinematic signature in K3-50A. However, this 
"fiow" is not aligned with the large scale ionized out fiow reported in earlier radio continuum 
and recombination line studies (lO's of arcseconds in extent). The NIFS kinematic signature 
is not symmetric about any detected point source in the NIFS image either. It is possible 
the bipolar kinematics arise from a lower mass protostar associated with the central star of 
K3-50A, and the associated material is ionized by the hot central source. A second broad 
component of Br7 emission is seen within a few tenths of an arcsecond of the continuum 
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peak. This may be hot gas escaping from the cavity resulting from the action of the ionizing 
radiation on the natal molecular cloud or large scale torus. 
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Fig. 1. — fC— band spectrum extracted from the central 2.5" diameter of the NIFS FOV. 
Differences between the object and sky frames resulted in a slight over-subtraction. Line 
identifications and strengths relative to Br7 are given in Table 1, and a number of the 
brighter lines are indicated in the figure. The H I Pfund series is seen near 23000 A and 
redward. 
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Fig. 2. — Continuum maps. Left panel: Log of the intensity grayscale for the 14.7 GHz radio 
continuum image from De Free et al. (1994, see their Figure 1). The angular resolution 
is 2.0" X 1.6" (beam FA=79°). The contours correspond to the NIFS continuum {right 
panel) where the peak emission in both images is made to match. Right panel: i^'— band 
continuum at 21700 A (19.2 A, 9 pixels wide) shown in grayscale (log intensity) and with 
the same contours plotted in the left panel. The field center is located at approximately 
RA(2000)=20/i 01m 45.71s, Dec(2000)= +33d 32' 42.1" (see text for a discussion of source 
position versus wavelength). The position angle on the sky was 70° (east of North); this image 
has been rotated so that North is up, and East is left. The red filled circles are sources from 
Hofmann et al. (2004) and the yellow diamonds are the mid-infrared sources of Okamoto et 
al. (2003); see text. The continuum peak appears as an unresolved point source near 0.2, 
0.6 and is assumed to be Source 7^1 of Hofmann et al. (2004) and 0KYM3 of Okamoto et 
al. (2003). The long line starting at Source #1 indicates the approximate direction of the 
ionizied outflow described by De Free et al. (1994). The second line is our estimate of the 
rotation axis of the HCO"*" torus described by Howard et al. (1997); see text. The torus is 
centered a few arcseconds N-W of the continuum peak. 
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Fig. 3. — The continuum subtracted Br7 map. The hues and points are the same as for 
Figure 2. Several of the continuum sources detected in Figure 2 are among the brightest 
Br7 sources. The brightest peaks appear angularly resolved, and it is not clear if there are 
associated point sources embedded within or if these are ionized clumps of higher density 
material. There are 14 contours overlaid on the image. The lowest contour is 1.48x10 
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cm~^ s~^ arcsec"^ and the contour spacing is the same. The peak is located at —0.6, 
and has a value of approximately 2.20x10 
is affected by He I blends; see text. 
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Fig. 4. — The ratio of He I to Br7 (grayscale and contours). The hnes and points are the 
same as for Figure 2. Left panel, 3p — 4s ^P°-'^S 21127 A line. The contours are spaced by 
0.01 and the peak ratio near Source ^1 is 0.055. The ratio (uncorrected for extinction) is 
uniform over most of the emission region and has a value of approximately 0.04 to 0.05. He I 
blends can affect the strength of Br7 as does a broad line component near the continuum 
peak; see text and also Table 2. The area near 0KYM4 has a smaller ratio, ~ 0.025, and this 
may be evidence of a lower excitation source. It also suggests the dominant continuum peak 
ionizes a volume corresponding to about 2" (see also Figure 10). Right panel. He I ^S-^P° 
20587 A line. In contrast to the 21127 A line, this ratio shows strong variations over the 
FOV: lower values (~ 0.25) at the continuum peak location and systematically increasing 
away from the peak to the south (~ 0.7) in the panel. 
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Fig. 5. — Higher ionization line maps. The hnes and points are the same as for Figure 2. 
[Fe III] 22184 A is shown on the left, and [Se IV] 22867 A is shown on the right. The 
ionization potential of Fe III is lower than that for Se IV showing that the excitation is 
higher closer to the continuum peak. The peak surface brightnesses are 6.6x10"^^ erg cm~^ 
s~^ arcsec"^ and 1.1x10"^^ erg cm~^ s~^ arcsec"^ for [Fe III] and [Se IV], respectively, and 
the contour intervals are evenly spaced at an interval 10 times less than the peak value. 



- 30 - 




1 -1 

Delta E (arcsec) 

Fig. 6. — Molecular hydrogen emission map in the v= 1-0 S(l) 21218 A line. The lines 
and points are the same as for Figure 2. The emission is clumpy and appears to trace the 
ionized gas (Figure 3). The peak emission is centered on the peak of the Br7 emission (shown 
as contours; see Figure 3). An equally strong peak is coincident with the continuum peak 
(Figure 2). The peak value of H2 emission is 1.7x10"^^ erg cm~^ s^^ arcsec"^ 
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Fig. 7. — Velocity (km s~^) for the principal component of the Br7 line in each NIFS spaxel 
derived from the data with a Gaussian fit to the line profile. The labeled lines and points 
are the same as for Figure 2. The mean Br7 velocity (Vlsr = —62.5 km s^^) of the principal 
component has been removed from the map. The bipolar structure does not line up with 
any known point source, and the axis does not align with the direction of the large scale 
ionized flow detected at radio wavelengths. The peak relative velocities are —24 km s~^ and 
+31 km s~^, but the color map is "stretched" slightly to show smaller velocities. 
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Fig. 8. — Detail of the line emission in the vicinity of Br7. The vertical hash marks identify 
He I lines (see Table 2). Five lines were fit simultaneously in each spectrum; see text. 
Left panel: emission from a 0.05" aperture located near Hofmann et al. (2004) source 7^6. 
Right panel: emission approximately 0.2" south and east from the continuum peak in an ~ 
0.05" aperture. The spectra have been continuum subtracted and scaled to the peak Br7 
flux. Each panel shows the observed spectrum in filled circles and the individual lines and 
continuum as solid black lines. The combined fit is shown as a solid red line. The broad line 
blended with Br7 in the right panel may be a distinct blue shifted Br7 component. It is 
much too strong to be due to He I under typical conditions; see text. The broad line in the 
right panel is about 120 km s^^FWHM compared to 70 km s^^ for the principal component 
of Br7 (i.e. the bright component). 
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Fig. 9. — Cloudy ionization model prediction of the He I 21127 A emission compared to Br7 
(color-scale and solid contours) for a model grid of T^q versus ionizing flux from the central 
star at the irradiated face of the H II region. The ratio increases toward higher Tes and 
ionizing flux. The contours have values of 0.01 to 0.08 and are spaced by 0.01. The solid 
red contour is the observed value (Table 2) increased by a few percent to 0.049 to account 
for differential extinction between the He I line and Br7. The dashed contours represent the 
ratio of Br7 to total dust emission and decrease toward lower T^s and higher ionizing flux. 
These contours range from 10~^ to 10""^, each separated by a factor of 10. The dashed red 
contour is the observed ratio of 1.8x10"^ (= 3.2x10"^^ erg cm~^ s~^ / 1.8xl0~^erg cm~^ 
s~^) with Br7 corrected for extinction; see text. 
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Fig. 10. — The ratio of He I to Br7 volume emissivities showing where the hnes are formed in 
the cloud. This plot is for a central object with T^s = 38000, LoglO(phi) = 14, and density 
10^ cm~^. Even though the 20587 A ratio changes significantly within the cloud, average 
values of the line ratio integrated through the cloud will be relatively large unless the density 
is low; see text and Figure 11. The 21127 A and 20587 A ratios integrated through the cloud 
are 0.052 and 0.601, respectively. The depth coordinate corresponds to the Cloudy model, 
converting cm to arcseconds assuming a distance of 7000 pc. 
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Fig. 11. — The predicted ratio of He I lines to Br7 from Cloudy for a range of nebular 
densities. These models correspond to a He abundance by number of 0.1 to H. The models 
were computed with a central star with T^s of 38000 K and the incident radiation field had 
a flux of ionizing photons, LoglO(phi) = 14.0. 



Table 1. Integrated Emission Line Ratios in K3-50A 



Line Identification^ 


Ao (A)b 


x{kY 


FWHM (km s"^) 


Ratio to Br7^ 


ri-2 v=i-U hy i) 


01010 


1 1 


o9 


U.UU/D ± U.UUUl 


[Fe ml ^H-'^G 

L J 


21457 


21454 


71 


0.0052 ± 0.0001 


Br7 (H I 7-4) 


21661 


21657 


71 


1.0000 ± 0.0002 


[Kr III] ^Pi-3P2 


21987 


21983 


57 


0.0041 ± 0.0001 


[Fe III] ^H-^G 


22184 


22183 


62 


0.0164 ± 0.0001 


H2 v=l-0 S(0) 


22233 


22230 


70 


0.0034 ± 0.0001 


[Fe III] =^H-3G 


22427 


22422 


54 


0.0086 ± 0.0001 


H2 v=2-l S(l) 


22477 


22473 


55 


0.0020 ± 0.0001 


[Se IV] 2P3/2-2Pi/2 


22867 


22863 


71 


0.0291 ± 0.0003 



'^Atomic emission line identifications were made using the on-line data base at NIST 
(Ralchenko et al. 2007) and Peter van Hoof's web page (http://www.pa.uky.edu/ pe- 
ter/atomic/) that uses energy levels from the NIST database. H2 lines are from Black 
& van Dishoeck (1987). 

^Rest wavelength, vacuum 

'^Observed wavelength, mean blue shift is —3.8 A, or approximately 1.8 pixels. 

"^Emission lines were extracted from a 2.5" diameter aperture near the center of field of 
view (see Figure 1). The Br7 line map is shown in Figure 3, but the line flux for this ratio 
was obtained by fitting multiple Gaussians to the He I blends and Br7; see text. 



Table 2. Integrated He I Lines in K3-50A 





Line Identification^ 


Ao (A)b 


A(Ar 


FWHM (km s"^) 


Ratio to Br7^ 


He I 4s 




20430 


20427 


79 


0.0055 ± 0.0001 


He I 2s 


- 2p ^S- P° 


20587 


20583 


67 


0.6778 ± 0.0002 


He I 4p 




20607 


20603 


89 


0.0067 ± 0.0002 


He I 3p 


- 4s - S 


21127 


21122 


70 


0.0462 ± 0.0001 


He I 3p 


- 4s ^P°-^S 


21138 


21134 


68 


0.0110 ± 0.0001 


He I 4p 


- 7s 3p°_3s 


21500 


21497 


70 


0012 ± 0001 


Hel 4d 


- 7/ 3D-3f° 


21614 


21609 


71 


0.0250 ± 0.0002 


Hel 4(i 


- 7/ iD-^F° 


21623 


21620 


70 


0.0100 ± 0.0002 


He I 4/ 


- 7g iF°-iG, 3f°-3G 


21647 


21642 


70 


0.0508 ± 0.0002 


He I 4/ 


- Id iF°-^D, 3f°-3D, Br7?'^ 


21655 


21650 


70 


0.0670 ± 0.0002 


Hel4rf 


- Ip 3D-3p° 


21821 


21817 


62 


0.0044 ± 0.0001 


He I 4p 


- 7rf ip°-iD 


21846 


21842 


62 


0.0023 ± 0.0001 


Hel 4p 


- 7s ip°-iS 


22291 


22287 


78 


0.0012 ± 0.0001 



^Atomic emission line identifications were made using the on-line data base at NIST 
(Ralchenko et al. 2007) and Peter van Hoof's web page (http://www.pa.uky.edu/ pe- 
ter/atomic/) that uses energy levels from the NIST database. Most of the fines have wave- 
lengths calculated from the term differences. 

^Rest wavelength, vacuum 

'^Observed wavelength, mean blue shift is —3.8 A, or approximately 1.8 pixels. 
^Emission lines were extracted from a 2.5" diameter aperture near the center of field of 



view (sec Figure 1). The Br7 line map is shown in Figure 3, but the hne flux for this ratio 
was obtained by fitting multiple Gaussians to the He I blends and Br7; see text. 

^Emission near the continuum peak is too strong to be due to He I; see text and Figure 8. 



